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Abstract 

A detailed analysis of ZZ production at the upgraded Fermilab Tevatron 
and the CERN Large Hadron Collider is presented for general ZZZ and ZZ^ 
couplings. Deviations from the Standard Model gauge theory structure for 
each of these can be parameterized in terms of two form factors which are 
severely restricted by unitarity at high energy. Achievable limits on these 
couplings are shown to be a dramatic improvement over the limits currently 
obtained by e^e~ experiments. 
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I. INTRODUCTION 



The Standard Model (SM) of electroweak interactions makes precise predictions for the 
couphngs between gauge bosons due to the non-abehan gauge symmetry of SU{2)l<S)U{1)y- 
These self-interactions are described by the triple gauge boson (trilinear) WWV, Z'yV, and 
ZZV (y = 7, Z) couplings and the quartic couplings. Vector boson pair production provides 
a sensitive ground for direct tests of the trilinear couplings. Deviations of the couplings from 
the expected values would indicate the presence of new physics beyond the SM. 

To date the SM has passed this rigorous test with no observed deviations from the SM 
values. The WWV couplings have been measured with an accuracy of 10 — 15% in W~^W~, 
single photon and single W production at LEP2 [||], and with 20 — 40% accuracy in Wj, WZ 
and W~^W~ production at the Fermilab Tevatron collider The Z'-fV couplings can be 

probed in Z7 production in e^e~ and in hadronic collisions. The LEP2 [0] and Tevatron 
experiments find the Z'-fV couplings to be smaller than 0.05 — 0.4, depending on the specific 
coupling considered. The ZZV couplings, on the other hand, are only loosely constrained 
at the moment through ZZ production at LEP2 Due to low event rates after branching 
ratios, or large backgrounds, ZZ production was not observed by the Tevatron experiments 
in Run I. 

In Run II of the Tevatron which will begin in 2001, an integrated luminosity of 2 — 15 fb~^ 
is envisioned 0, and a sufficient number of ZZ events should be available to commence a 
detailed investigation of the ZZV couplings. At the CERN Large Hadron Collider [(LHC), 
pp collisions at ^/s = 14 TeV 0], one can imagine that the measurement of these couplings 
reaches the 0.1% level of current precision electroweak data. In this paper we study the ca- 
pabilities of future hadron collider experiments to test the ZZV vertices via ZZ production. 
In the past, the reaction pp^ ^ ZZ for non-zero ZZV couplings has only been studied in 
the approximation where the Z bosons are considered as stable final state particles |lTO| , p 



We go a step further and include Z decays with full decay correlations, finite Z width effects 
and time-like virtual photon exchange in our analysis. 

Two ZZZ couplings, and two couplings, are allowed by electromagnetic gauge 

invariance and Lorentz invariance for on-shell Z bosons. We discuss the properties of 
these couplings in Section ||, where we also derive unitarity bounds for the form factors 
associated with the ZZV vertices. The SM is assumed to be valid in our calculations except 
for the ZZV anomalous couplings; Vff couplings and strong interactions of SM particles 
remain unchanged. 

Our analysis examines the observable final state signatures at hadron colliders, ZZ — >■ 
diiit^^i 1 h, £2 = e, /i) and z/i/jj. In Section |TTl| we provide details 

of the signal and various backgrounds and discuss the signatures of anomalous ZZZ and 
ZZ'y couplings. Besides the ZZ invariant mass distribution and the Z boson transverse 
momentum distributions, the azimuthal angle between the Z boson decay fermions, A$, and 
their separation in the pseudo-rapidity - azimuthal angle plane, AR, are sensitive indicators 
of anomalous couplings. The A$ distribution may be useful in discriminating different types 
of ZZV couplings, and in determining their sign. In Section |^ we derive sensitivity limits 
for anomalous ZZV couplings for various integrated luminosities at the Tevatron and LHC 
and discuss the results. Finally, in Section [V| we compare the expectations for Tevatron 
Run II and the LHC with the current LEP2 limits and expectations for an e~^e~ Linear 
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Collider. In Section we also present our conclusions. 



II. ZZZ AND ZZ7 ANOMALOUS COUPLINGS 

In the SM, at the parton level, the reaction — ZZ proceeds by the Feynman diagrams 
of Fig. m. Including the anomalous couplings under discussion requires the addition of the 
graphs of Fig. 0. In the massless fermion limit, a reasonable approximation for hadron 
collider processes, the most general form of the Z^{q2) V^{P) (V = Z, 7) vertex 

function (see Fig. |^ for notation) for on-shell Z's which respects Lorentz invariance and 



electromagnetic gauge invariance may be written as |T2 

P2 _ Ml 



zzv - ^ T72 



Ml 



(1) 



where Mz is the Z-boson mass and e is the proton charge. The overall factor (P^ — My) in 
Eq. ([l|) is a consequence of Bose symmetry for ZZZ couplings, while it is due to electromag- 
netic gauge invariance for the ZZ7 couplings. The couplings fY {i = 4, 5) are dimensionless 
complex functions of g| and P^. All couplings are C odd; CP invariance forbids fY 
and parity conservation requires that vanishes. Because and fj are CP-odd, contri- 
butions to the helicity amplitudes proportional to these couplings will not interfere with the 
SM terms. In the static limit, f2 corresponds to the anapole moment of the Z boson |]T3 



In the SM, at tree level, fj = fY = 0. At the one-loop level, only the CP conserving 



couplings fY receive contributions. Numerically, these contributions are of O(10~^) 
Loop contributions from supersymmetric particles and additional heavy fermions produce 
ZZV couplings of a similar magnitude [jl4|. If the Z bosons are allowed to be off-shell, five 



additional ZZZ couplings, and five additional ZZ'~f couplings contribute |]T5|. For these 
couplings, the factor (P^ — My) in the vertex function is replaced by {q\ — q^)- The effect 
of these couplings thus is strongly suppressed and we shall ignore them in our discussion. 

It should be noted that the two ZZ'-j couplings contributing to ZZ production are 
completely independent of the four Z7Z couplings which appear in Z7 production (assuming 
that the Z-boson is on-shell). If all three vector bosons in the vertex function are off-shell, 
there are seven couplings altogether. Four of them survive in Z7 production, and two in 

ff-^zz. 

The parton level di-boson production cross sections with non-SM couplings manifestly 
grow with the parton center of mass energy V^- S'-matrix unitarity restricts the ZZV 



couplings uniquely to their SM values at asymptotically high energies [T^. This requires 
that the couplings fY possess a momentum dependence which ensures that the fYiQi, P"^) 
vanish for any momenta much larger than Mz- For ZZ production, g^, g| ~ M| even 
considering finite Z width effects, but P"^ = s may be quite large at the hadron colliders 
under consideration. In order to avoid unphysical results that would violate unitarity, the 
s dependence thus has to be taken into account. To parameterize the s dependence of the 



form factor, we will use a generalized dipole form factor |T7 

fV 

[l + s/Al^)- 



/^(^1 = 7^^^ (^ = 4,5), (2) 
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where ApF the form factor scale which is related to the scale of the new physics which is 
generating the anomalous ZZV couplings. 

The values = f^{M^, M|, 0) of the form factors at low energy (s = 0) and the power 
of the form factor, n, are constrained by partial wave unitarity of the inelastic ZZ produc- 
tion amplitude in fermion antifermion annihilation at arbitrary center-of-mass energies. In 
deriving unitarity limits for the /j^'s, we follow the strategy employed in Ref. |jl8|. The 
anomalous contribution to the 

f{cr)f{a) ^ Z{X,)Z{X2) (3) 
helicity amplitudes may be written as 



(4) 



where are the conventional d-functions /5 = (1 — 4M|/s)^/^, a, a and Ai,A2 are 
the helicities of the incoming fermion pair and outgoing Z pair, respectively. g^J^^^ is the 
coupling of the s-channel vector boson to the incoming fermion pair, B is the center of mass 
scattering angle, and the A^'s are the reduced amplitudes given by (see also Ref. |2( 



A 



o± 



2Mz 

s 



2M 



^00 



4^ 



0. 



Examination of the J = 1 partial wave amplitude produces the desired unitarity bounds. 



A1A2 



1/2 



< 



1 



a/33/2 



- ( 3 — 6 sin^ 9w + 8 sin"^ 9w 



1/2 Ml 



A2I 



A1A2 



1/2 



< 



3 . . . M| 



a/53/2 y 10 



— sin 9w cos 9 



w 



(5) 

(6) 
(7) 

(8) 
(9) 



where 9\y is the weak mixing angle and a the QED fine structure constant. 

By substituting Eq. (H) and assuming that only one coupling is nonzero at a time, we 
find the following unitarity bounds for KpF 3> M^: 



1/40,50 1 



1 

< - 

a 



^ ( 3 — 6 sin^ 6'vi/ + 8 sin^ 9w 



V2 fMz\^ 
Aff) 



(n-3/2) 



1/ 



40,50 1 



< 



a V 10 



3 . , n fMz 
sm 9w cos 9w 



A 



FF 



■ n 



jn-l 



{n-3/2) 



(10) 



(11) 



Tree level unitarity is satisfied throughout the entire s range when these limits are observed. 
For the more likely case that several anomalous couplings contribute, cancellations may 
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occur and the bounds are weaker than those hsted in Eqs. (|T^) and (|Tl|). From the n 
dependent factors in ([ToD and (11) one observes that n > 3/2 is necessary in order to 



satisfy unitarity. This is a direct consequence of the high energy behavior of the anomalous 
contributions to the ZZ hehcity amphtudes, which grow hke {y/S/Mzf. In the following 
we shall assume that n = 3. Selecting an exponent sufficiently above the minimum value 
of 3/2 ensures that the ZZ differential cross section stays well below the unitarity limit 
at energies a/5 ^ Aff ^ Mz, where novel phenomena such as resonance production, or 
multiple weak boson production, are expected to dominate. For the form factor scale we 
choose Aff = 750 GeV at the Tevatron and Aff = TeV at the LHC in our numerical 
simulations. 



III. SIGNATURES OF ANOMALOUS ZZV COUPLINGS 

In this section we discuss characteristics of the signal and backgrounds of anomalous 
ZZV couplings. For simplicity, we consider only real ZZV couplings. We consider four 
signatures of ZZ production: decays to four leptons; two leptons and missing energy; two 
leptons and two jets; and two jets plus missing energy. Due to the overwhelming four 
jet QCD background, decays where both Z bosons decay hadronically are not considered 
here. We calculate the SM signal, the signal with anomalous couplings, and the significant 
backgrounds via full tree level matrix elements for the subprocess in question, each of which 
is discussed in detail below. 



A. General considerations 

Our calculation is carried out at the tree level. We compute the qq ZZ —>■ 4 fermion 
helicity amplitudes in the double pole approximation which ignores contributions from non- 
resonant diagrams except for contributions from time-like photon exchange diagrams, using 
the method described in Ref. Decay correlations, finite Z width effects and contri- 



butions from time-like photon exchange are taken into account. Cross sections and dynam- 
ical distributions are evaluated using a parton level Monte Carlo program. 

To simulate the effect of next-to-leading-log (NLL) QCD corrections we multiply the 
differential cross section with a simple if-factor which depends on the final state considered. 
A more detailed discussion of how NLL QCD corrections affect the four different final states 
is presented in the following sections, where we discuss each final state in turn. Gluon fusion, 
gg ZZ, contributes about 1% (15%) to the cross section at the Tevatron (LHC) [^]. We 
do not include the contribution from gluon fusion into our analysis. 

To examine the effects of anomalous couplings on observables we simulate pp {pp) colli- 
sions at ^/s = 2 TeV {^/s = 14 TeV) for the Tevatron (LHC). For all our numerical results 
we have chosen the set of SM input parameters to be: sin^ 9w = 0.2310, Mz = 91.187 GeV, 



and a{Mz) = 1/128.93 ^3]. For all processes which depend on the QCD coupling constant. 



we choose the value of the strong coupling constant to be as{Mz) = 0.118. We employ 



CTEQ4L parton distribution functions ||2J] for all calculations, selecting the value of the 
factorization scale to be /i/ = Mz- 
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As finite detector resolution can have a sizable effect on cross sections and thus the 
number of events accepted into the data set, to make our calculations realistic we must take 
into account some minimal detector response. We accomplish this via Gaussian smearing of 
the four momenta of the outgoing particles according to detector expectations. For Gaussian 



smearing at the Tevatron we use the expected values of the upgraded CDF detector |25 
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E 

For the LHC we take the expected values for the ATLAS detector |2 
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E 
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Here, rj is the pseudo-rapidity, E {Et) is the energy (transverse energy) measured in GeV, 
and the © sign symbolizes that the two terms are added in quadrature. 

In all cases, the missing momentum in an event is taken as the negative vector sum of the 
smeared four momenta of all observable final state particles. This does ignore the effects of 
additional soft activity that will affect this distribution in experiment, but for our purposes 
here may be safely neglected. 

The geometric and kinematic acceptance of detectors, i.e. the ability to observe and 
properly identify final states particles, is simulated in our calculations by cuts imposed on 
observable particles in the final state. At the Tevatron (LHC) we require {i = e, /i): 



Tevatron 



LHC 



PT{i) > 15 GeV 

\vm < 2.5 
PtU) > 20 GeV 
\vij)\ < 2.5 
AR{ej) > 0.6 
AR{jj) > 0.6 
6j, > 20 GeV for ZZ 
^j, > 60 GeV for ZZ - 



vv 



PT{i) > 15 GeV 

\rj{£)\ < 2.5 
PtU) > 30 GeV 

\v{j)\ < 3 
AR{ij) > 0.6 
ARijj) > 0.6 
:^r^ > 50 GeV for ZZ . 
^j, > 60 GeV for ZZ - 



where pt is the transverse momentum and 



AR 



{A<^f + {At^Y 



1/2 



(12) 
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is the separation in the pseudorapidity - azimuthal angle plane. is the missing transverse 
momentum resulting from the nonobservation of the neutrino pair. 

In addition, we require the C.^^~ and two jet invariant masses to be within ±15 GeV of 
the Z boson mass [p^ , ^ : 



76 GeV < m(^+r)< 106 GeV, (13) 
76 GeV < m(jj) < 106 GeV. (14) 

These cuts help suppress contributions from non-resonant Feynman diagrams and, in the 
jet case, the background from Z + 2 jet production. Additional cuts which are imposed to 
reduce backgrounds for individual final states will be discussed when we consider the specific 
final state to which they apply. 



B. The qqqq channel 

The first ZZ decay channel we consider, ZZ iiiii2i2 (^i? ^2 = e,/i) is observation- 
ally the cleanest as it is essentially background- free. However, it does suffer from a small 
event rate due to a tiny branching ratio of B{ZZ ifi^i2^2 ) = 0.0045 if both electron and 
muon final states are considered. In the following we concentrate on the ZZ e^e~ ix^ ijT 
channel. Results for the e^e~e^e~ and final states can be obtained by dividing 

all cross section results by two. However, for these final states one has to take into account 
the combinatorial background originating from not being able experimentally to distinguish 
identical charged leptons. All results presented in this section include a /^-factor of 1.28 



(1.34) at Tevatron (LHC) energies to approximate the effect of NLL QCD corrections [29 
Inclusive NLL QCD corrections to ZZ production are known to modify the shape of distri- 
butions only insignificantly. As mentioned in Sec. [HI A| , we do not include the non-resonant 



Feynman diagrams which contribute to qq — e~^e~ fi~^n~ in our calculation. Requiring the 
invariant mass of the lepton pairs to be in the vicinity of the Z mass, (see Eq. (p!3D), these 
diagrams contribute less than a few percent to the differential cross section. Imposing the 
cuts listed in Sec. |III A| , one obtains a ZZ — >■ £~^£~£~^i~ (£ = e, /i) cross section of 3.85 fb 



(22.5 fb) at the Tevatron (LHC). For 2 fb~^ in Run II, only a few ZZ 4 lepton events are 
therefore expected within the framework of the SM. 

Similar to the WWV and Z'jV couplings [[I2| , p!7| , the effects of anomalous ZZV cou- 
plings are enhanced at large energies. A typical signal of nonstandard ZZZ and ZZ'j 
couplings thus will be a broad increase in the ZZ invariant mass distribution, the Z trans- 
verse momentum distribution and the pt distribution of the Z decay products. The Pt{Z) 
and prifJ') distributions for the Tevatron {pp collisions at = 2 TeV) and the LHC are 
shown in Figs. § and ^, respectively. Results are shown for the SM and two ZZV couplings, 
f^Q = 0.3 (0.02) and f^Q = —0.3 (—0.02) at the Tevatron (LHC). Here, and in all subsequent 
figures, only one ZZV coupling is allowed to be non-zero at a time. 

Terms proportional to f]^ and in the matrix elements have identical high energy be- 
havior. Differences in the differential cross sections at high energies between ZZZ and ZZ7 
couplings are thus controlled by the Zff and 7// couplings, and by the parton distribution 
functions. At the Tevatron these result in differential cross sections which differ by only a 
few percent for s ^ M| if \f^\ = |/jo| {i = 4,5). Slightly larger differences are observed at 
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intermediate energies and transverse momenta. Since fY violate CP conservation, terms in 
the helicity amplitudes proportional to those couplings do not interfere with the SM terms. 
Cross sections thus are independent of the sign of fY- Interference effects between the 
anomalous and SM contributions to the helicity amplitudes, however, do occur for . The 
magnitude of the interference effects in the transverse momentum and ZZ invariant mass 
distributions unfortunately is small. 

While it would be difficult to discriminate between the various ZZV couplings in the 
^zz or pt distributions, it should be easy to distinguish between the SM Higgs boson which 
decays in a pair of Z bosons, H —>■ ZZ, and anomalous ZZV couplings. This will be impor- 
tant at the LHC, where final states resulting from pp ZZ are prime Higgs boson search 
channels. Anomalous gauge boson couplings lead to a broad increase in the differential cross 
section at large energies and transverse momenta, whereas a scalar resonance produces a 
Breit-Wigner resonance in the mzz distribution and a Jacobian peak in the Pt{Z) spectrum. 
In addition, the correlation of the angular distributions of the Z decay leptons in the Z bo- 
son rest frames may be used to discriminate between a Higgs boson and ZZV couplings. Z 
bosons originating from Higgs boson decay are mostly longitudinally polarized |30|, whereas 
anomalous ZZV couplings lead to one transversely polarized and one longitudinally polar- 
ized Z boson (see Eqs. (^) and (|^)). Since the Z boson coupling to charged leptons is almost 
purely axial vector like, transversely polarized Z bosons lead to a angular distribution for 
the Z decay leptons which is proportional to (1 + cos^ 9*), where 6* is the polar angle in the 
Z boson rest frame with respect to the flight direction of the Z boson in the ZZ rest frame. 
On the other hand, the angular distribution for longitudinal Z's is proportional to sin^ 6* . 

In order to distinguish and and to determine the sign of the Ai? distribution 
and the distribution of the opening angle in the transverse plane, A$, of the t^^~ pair origi- 
nating from the decay of a Z boson may be helpful, if deviations from the SM predictions are 
found in the mzz and the transverse momentum differential cross sections. Figure ^ shows 
the Ai?(/i"'"/i~) and A$(/i"'"/i~) distributions for ZZ e+e~/i+/i~ at the Tevatron in the 
SM and for non-standard ZZZ couplings. Similar results are obtained for the correspond- 
ing distributions of the e^e~ pair, and for the ZZ7 couplings fj^- The SM AR{fi'^ fi~) and 
A^{fi^ differential cross sections are dominated by the threshold region, ^ 2Mz, 
where the Z boson momenta are small and the decay leptons tend to be back-to-back, i.e. 
the distributions are strongly peaked at AR ^ 3 and A$ = 180°. Anomalous couplings 
affect the cross section mostly at large Z-boson transverse momentum. Due to the Lorentz 
boost, the relative opening angle between the leptons originating from Z decreases 
with increasing pt- The deviations due to non-standard ZZV couplings in the AR{n'^ n~) 
and A$(yU+;U^) distributions are therefore concentrated at rather small values. Figure ^ 
demonstrates that the A$(/i+/i~) distribution would be particularly useful in separating 
the individual ZZV couplings. The shapes of the A$(/i"'"/i~) distributions for = 0.3, 
f§) = 0.3 and = —0.3 differ considerably, and for a sufficient number of events it should 
be possible to distinguish fY from f^, and to determine the sign of f^- Unfortunately, both 
the AR{fi^ fi^) and the A$(yU~^/i~) differential cross sections are useless in distinguishing ff 
from /7 {i = 4,5). As for the ZZ invariant mass distribution and the transverse momen- 
tum distributions, the differential cross sections for ff = differ very little in shape and 
magnitude. 

At the Tevatron, the small number of ZZ — > 4 leptons events will limit the usefulness 
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of the Ai?(£"'"£~) and A^{£~^£~) distributions. At the LHC, the expected number of events 
is much larger; however, the magnitude of the interference effects between the SM and 
the anomalous contributions to the helicity amplitudes is significantly smaller than at the 
Tevatron. This can be easily understood from the high energy behavior of the anomalous 
contributions to the differential cross section. The differential cross section is proportional 
to the squared amplitude, which contains the SM terms, terms linear in the anomalous 
couplings and terms which are quadratic in the ZZV couplings. The terms linear in the 
anomalous couplings originate from the interference between the SM amplitude and the 
anomalous contributions, and are proportional to (-^/I/M^)^. Terms quadratic in the ZZV 
couplings on the other hand are proportional to (s/M|)^ and thus grow much faster with s. 
Due to the much higher parton center of mass energies accessible at the LHC, interference 
effects thus play a smaller role than at the Tevatron. The AR{fi~^ fi~) and the A<I>(yLi"'"yLi~) 
distributions at the LHC for the SM, f^Q = 0.02, /g^ = 0.02, and fi, = -0.02 are shown 
in Fig. 1^. Both distributions are very insensitive to the sign of The shape of the 
A<l>(/i+/i") distribution differs for ff^ and f§) for A<l>(/i+/i") < 20°, whereas the Ai?(/i+/i") 
distributions for and values of equal magnitude are almost identical. It would thus 
be challenging to discriminate between and f^, and to determine the sign of f^, using 
the Ai?(£+£-) and the A$(£+£-) distributions in the ZZ ^ 4 leptons mode at the LHC. 
Since events in regions where anomalous ZZV couplings have a significant effect originate 
from higher values of VI, the typical separation and the opening angle in the transverse 
plane are significantly smaller than at the Tevatron. 



C. The i+ruu channel 

In contrast to the ZZ — ^ 4 leptons mode which is almost background free, there are 
several potentially important background processes if one of the two Z bosons decays into 
neutrinos. The advantage of the ZZ — >■ i^i~ui> channel, observable as i^i~^rp, is its larger 
branching fraction. Summing over the three neutrino species, the raw number of ZZ — >■ 
signal events before cuts are implemented is about a factor 6 larger than the number 
of ZZ 4 leptons events. 

Calculation of the signal is similar to that for the four lepton channel, with the change 
of one Z coupling from gzee to gzuu- There is no photon interference for the Z that decays 
to neutrinos. Since it is not possible to restrict the invariant mass of the neutrino pair to 
be in the vicinity of the Z resonance, one has to verify that the non-resonant Feynman 
diagrams ignored in our calculation do not significantly change the differential cross section. 
There are 10 Feynman diagrams contributing to qq — ^ £^£~ueii>ei, i ^ and 19 contributing 
to qq i'^i'h'ei'i in the SM. The graphs of the qq W^W~ (i^(i~UfPt background. 



which is described in more detail below, are contained in this set. Using madgraph ||31 



and the helas |^ library, we have calculated the SM cross section for pp e^e'vu 
including the full set of contributing Feynman diagrams and summing over the three neutrino 
species. In Fig. ^we compare the transverse momentum distribution of the e^e~ system for 
pp —>■ e^e^vv at the Tevatron resulting from the full set of tree level Feynman diagrams with 
the distribution obtained using the subset of diagrams contributing to qq —>■ ZZ e^e^vu 
and qq W^W~ — > e+e~z/eZ^e in the double pole approximation. Figure |a displays the 
individual pt distributions, whereas part b) of the figure shows the ratio of the differential 
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cross sections. In addition to the cuts listed in Sec. pi7\| we have imposed a > 20 GeV 
cut in Fig. |^. The non-resonant diagrams are seen to reduce (enhance) the rate by about 
5% for pT{e~^e~) < 80 GeV (pT(e'''e~) > 80 GeV). This is significantly smaller than other 
theoretical uncertainties such as the factorization scale uncertainty in our calculation. In the 
following we shall therefore ignore effects of the non-resonant diagrams in channels where 
one of the Z bosons decays into neutrinos. 

The most important background processes to ZZ — > production are tt — > 

W^W'bb, standard electroweak W^W~ + X production with W^W~ and 
Z(— 5> + 1 jet production with the jet rapidity outside the range covered by the de- 

tector, thus faking missing px- We will call this last process the "Z + 1 jet" background. 
We have calculated the ti background using standard helicity amplitude techniques, fully 
including the subsequent decays t — > Wb and W ^ iu and all decay correlations. Finite 
width effects for the top quarks and W^s are included. Jets (partons) with AR{jj) < 0.4 
are merged into a single jet. We do not decay the bottom quarks explicitly, but do include a 
parameterized energy loss distribution to make a more realistic simulation of observed final 
state momenta and overall missing momentum. For W^W~ +X production we make use of 
the calculation described in Ref. |Q. For a realistic assessment of the Z + 1 jet background, 
a fuU-fiedged Monte Carlo simulation is required. Here, for a first rough estimate, we use a 
simple parton level calculation. For a jet, i.e., a quark or gluon, to be misidentified as j^^, at 
the Tevatron (LHC), we shall require that the jet pseudorapidity be \r]{j)\ > 3.5 (|r7(j)| > 5). 

Additional backgrounds originate from bb production and Z t^t^ decays. These 
backgrounds can be suppressed to a negligible level by requiring the angle in the transverse 
plane between a charged lepton and the missing transverse momentum to be between 20° 
and 160°, if < 50 GeV. 

Subsequently, in this section we shall focus on the e^e~^j, final state. Virtually identical 
results are obtained for ZZ fi^fi^^j,. For reasons which will become clear shortly, we do 
not include a /^-factor for the signal cross sections in the figures shown in this section. In 
addition to the cuts specified in Sec. pi A| , we impose a i^^, > 20 GeV {-^rp > 50 GeV) cut 
at the Tevatron (LHC), and the cut on the angle in the transverse plane between a charged 
lepton and the missing transverse momentum discussed above. 

In Fig. 1^ we show the transverse momentum distribution of the e~^e~ pair for ZZ 
e~^e~'^rp in the SM (solid curve) and for = 0.3 (dashed line) at the Tevatron, together 
with the differential cross sections of the ti (dotted line), W^W~ (dot-dashed line) and 
Z + 1 jet (long-dashed line) backgrounds. The cut on the angle in the transverse plane 
between a charged lepton and the missing transverse momentum is responsible for the slight 
dip in the ZZ differential cross section curves at pxic^e^) ~ 50 GeV in Fig. 

One observes that all backgrounds shown in Fig. |^ are significantly larger than the 
ZZ signal at small values of pT(e^e~). Because of kinematical constraints, however, the 
Z + 1 jet background drops very rapidly with p^- In a more complete treatment in which 
soft gluon and/or quark radiation and hadronization effects are included, one expects that 
the pT{e~^e~) distribution will be somewhat harder for the Z + 1 jet background, especially at 
high transverse momentum. The Z -|- 1 jet background sensitively depends on the rapidity 
cut above which a jet is assumed to be misidentified as Our assumption that jets 
with |?7(j)| > 3.5 will fake missing transverse momentum is probably conservative, and the 
Z + 1 jet background may well be significantly lower than shown in Fig. |^. 
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The W^W~ background exceeds the ZZ signal cross section for pT(e"'"e~) < 80 GeV. 
Since the tail of the pr(e"'"e~) distribution for W^W~ production is very sensitive to NLL 
QCD corrections, we show the NLL pT{e~^e~) distribution in Fig. The differential cross 
section of the tt background is larger than the SM signal for e^e~ transverse momenta as 
large as 200 GeV, and may thus reduce the sensitivity to anomalous ZZV couplings. We 
do not show the distributions for the bb or Z ^ t^t^ backgrounds, as they are negligible 
after the aforementioned angular cut. 

The ti background can be virtually eliminated by requiring that no jets with pxii) > 
20 GeV and \ri{j)\ < 3.5 are present. Such a jet veto also reduces the W~^W~ +X background 
at large transverse momenta. This is shown in Fig. ^d. The remaining W'^W" background 
will only marginally affect the sensitivity to ZZV couplings. The jet veto also significantly 
reduces the size of the NLL QCD corrections for the ZZ signal, justifying our procedure of 
not including a i^-factor in the signal cross section in this section. 

The pT{e^e~) distributions for signal and background processes at the LHC, requiring 



that no jets with prij) > 50 GeV and \ri{j)\ < 5 are present, are shown in Fig. |T0[ Since 
the ti production cross section is more than two orders of magnitude larger than at the 
Tevatron, the tt background is significant even when a jet veto is required. For the cuts 
imposed, the tt background (dotted line) exceeds the signal for pj'{e~^e~) < 140 GeV but 
is negligible for large transverse momenta. Both the W~^W~ and tt backgrounds could 
be reduced somewhat by choosing a smaller e~^e~ invariant mass window. Due to the more 
severe missing transverse momentum cut and the improved rapidity coverage of the hadronic 
calorimeters of the LHC experiments, the Z + 1 jet background is very small. 

As in the ZZ 4 leptons channel, the AR{e~^e~) and A$(e"*"e~) distributions would 
be useful at the Tevatron in distinguishing between and , and in determining the 



sign of f^. As shown in Fig. ITI|, the W^W background is negligible in those regions 



of Ai?(e^e^) and A$(e^e^) where the contributions from anomalous ZZV couplings are 
most pronounced. The jet veto imposed in Fig. |TT| renders the ti background negligible (see 
Fig. Pd). The '^j, cut removes events where the Z bosons are produced right at threshold 
and thus causes the peak in the AR{e^e~) (A$(e"'"e~)) distribution to shift from ^ 3 (180°) 
to ~ 2.6 (~ 140°). It has a similar effect on the W^W~ + jet background. Note that the 
W pair production background vanishes for Ai?(e"'"e^) < 1.4 and A$(e^e^) < 90°. The 
Z + 1 jet background is not shown in Fig. |TT]to avoid overburdening the figure. Qualitatively, 
the Ai?(e"'"e~) and A$(e"'"e~) distributions of the Z + 1 jet background are similar to those 
of the W^W~ + jet background. 

For completeness, we show the Ai?(e"'"e~) and A$(e"'"e~) distributions at the LHC in 
Fig. |T2|, imposing the same jet veto requirements as in Fig. Due to the higher missing 
transverse momentum cut, the peak in the SM ZZ AR{e~^e~) (A$(e+e~)) distribution is 
shifted to AR{e+e-) ^ 2 (A$(e+e-) 100°). The AR{e+e-) distributions of the ti and 
W~^W~ + jet background peak at similar values. Both backgrounds are negligible for 
AR{e~^e~) < 1.4. In the SM, the dominant helicity at high energies in uu — ^ W^W~ 
{dd ly+iy-) is Xw± = =Fl (Avi/± = ±1) il- Because of the V - A nature of the Wiu 



coupling, the charged leptons in W^W~ — ^ e+e~z/et'e tend to be emitted either both into {dd 
annihilation), or both against the flight direction of their parent W boson {uu annihilation). 
If a jet veto is imposed, the and W~ in W pair production are almost back to back in 
the transverse plane, and so are the W decay leptons. The A$(e"^e~) distribution of the 
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W^W~ + jet background thus peaks at a significantly larger angle than that of the ZZ 
signal and of the tt background (see Fig. |T2|b). While the tt background is smaller than the 
SM ZZ signal for A$(e+e~) < 50°, it is not negligible in this region. 



D. The semi-hadronic channels 

The decay modes where one of the two Z bosons decays hadronically have much larger 
branching fractions than the ZZ 4 leptons and the ZZ l^l^vv channels, but also 
much higher backgrounds. Nevertheless, these channels may be useful in searching for ZZV 
couplings: both CDF and D0 [|] have successfully used the WW, WZ — > £z/jj channels 
to constrain anomalous WWV couplings in the past. 

The main background sources are QCD Z+2 jet {''Zjj") production and WZ production 
where the W decays hadronically into a pair of jets. In both cases the jet pair invariant mass 
is constrained to be near the Z pole, Eq. (|Tip. The range of jj invariant masses considered 
here roughly corresponds to Mz ± 2ajj where ajj is the two jet invariant mass resolution 
of the detector, which typically is ajj = 5 — 9 GeV p6i^ . The Z boson decays to either 
e~^e~ , fi~^fi~ or invisibly to neutrinos. 

The Zjj background consists of QCD real-emission corrections to Z production. These 
subpro cesses include p6| 



qg^qgi+r, qq'^qq'tr, (15) 

which are dominated by t-channel gluon exchange, and all crossing-related processes, such 
as 

qq -> ggtr, gg qqi^i~ . (16) 

Similar to the treatment of the signal processes, we use a parton-level Monte-Carlo program 
based on the work of Ref. to model the QCD Zjj background. All interference effects 



between virtual photon and Z-exchange are included for charged lepton final states; for 
final state neutrino pairs there is no photon contribution, as running at one loop order is 
included, correcting the contribution from each phase space point from the input value of 
as{Mz) = 0.118. To compute the WZ background we make use of the calculation presented 
in Ref. 0. 

Other potentially dangerous background sources are tt and Wjj production. The Wjj 
background contributes only to the i^i'jj final state. For ZZ —>■ i'^i~jj, the ti —>■ uii' iyibb 
background can be reduced by requiring that the missing transverse momentum is ^j, < 
20 GeV at the Tevatron, and '^j, < 40 GeV at the LHC. The higher value for the LHC is 
motivated by pile-up effects due to the large number of interactions per beam crossing at 
the LHC design luminosity of £ = 10^^ cm^^ s~^. Pile-up effects amplify small momentum 
imbalance effects due to mismeasurements in the underlying event structure and can lead to 
a significant fake missing transverse momentum. Backgrounds from W^W~jj production 
are significantly smaller than those from tt production after imposing a '^j, veto, and are 
ignored in the following. 

For ZZ i^vjj, suppression of the ti —>■ i^v^i'Vfbh and W{—^ ^^)jj backgrounds is 
possible by requiring that there are no charged leptons withp7^(£) > 10 GeV and \rj{t)\ < 2.5 
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present in the event. The tt iujjbb background can be further reduced by imposing the 
constraint that the event does not contain more than two jets satisfying the pt and pseudo- 
rapidity cuts described in Sec. [ill A| . The Wjj background is calculated using the results of 
Ref. ifl. 

In the following, for ZZ ^^^^jj, we shall sum over electrons and muons in the final 
state and impose a j^^, < 20 GeV {^j, < 40 GeV) cut at the Tevatron (LHC) in addition to 
the charged lepton and jet pt, pseudorapidity and invariant mass cuts specified in Sec. [ill A . 
For ZZ —y jyjyjj, we require a missing transverse momentum of > 60 GeV and no charged 
leptons with Pt(^) > 10 GeV and < 2.5. In addition, the number of jets which satisfy 

the cuts detailed in Sec. [Ill A[ , rij, has to be Uj = 2. The rather high cut ensures that 
backgrounds from heavy quark production and three jet production, where the rapidity of 
one of the jets is outside the range covered by the detector, are sufficiently suppressed. At 
the LHC, a > 60 GeV cut may well be too low to trigger on jj;^rp events, especially at 
high luminosities. However, since anomalous ZZV couplings lead to deviations only at large 
transverse momenta, raising the ^j, cut to 100 GeV or even 200 GeV at the LHC will have 
very little impact on the sensitivity of the ZZ vvjj mode to ZZV couplings. Finally, 
all signal and background cross sections presented in this section are calculated at leading 
order. 

The pt{^^^~) distribution for ZZ — > i^i^jj is shown in Fig. 13a for the Tevatron, while 
Fig. |l^b illustrates the Tevatron prUj) distribution for uujj events. We display the SM 
cross section (solid line) together with the main backgrounds. We also show the ZZ cross 
section for = 0.3 (dashed line). The "kink" in the WZ and ZZ differential cross sections 
at Pt ~ 250 GeV is due to the AR{jj) > 0.6 cut which becomes effective only at sufficiently 
high transverse momenta. Because of the missing transverse momentum (charged lepton) 
veto for ZZ i^^^'jj {ZZ vvjj), the ti background is negligible at the Tevatron. The 
Wjj background, which contributes only to the vvjj final state, is considerably larger than 
the SM ZZ signal only for small values of prijj), thus it will not affect the sensitivity to 
anomalous couplings significantly. The WZ differential cross section (dotted line) is very 
similar to that of the SM signal over most of the px range considered, whereas the Zjj 
background dominates, overwhelming the signal even at very high transverse momenta. 
Its size is uniformly about one order of magnitude larger than the SM ZZ signal. It will 
therefore be very difficult to observe ZZ production in the semi-hadronic channels, if the SM 
prediction is correct. However, for sufficiently large anomalous ZZV couplings, the ZZ cross 
section exceeds the background at large transverse momenta. The semi-hadronic channels 
therefore may still be useful in obtaining limits on the ZZV couplings at the Tevatron. 

The case is much worse for the LHC, as shown in Fig. O. There the rate for Zjj events 
is almost two orders of magnitude greater than the SM signal, and approximately one order 
of magnitude greater than that for ZZ events for moderate values of [i = 4,5). The 
Wjj background is a factor 3 to 10 larger than the SM ZZ vvjj signal. Although the 
ti background is significantly larger than at the Tevatron, it has almost no effect on the 
sensitivity of the jj'^rp final state to ZZV couplings. 

It should be noted that NLL QCD corrections could worsen the signal to background 
ratio. QCD corrections enhance the ZZ signal cross section by about a factor 1.3 The 



full NLL QCD corrections to Zjj and Wjj production are currently not known. However, 
the QCD corrections to Zhh production in the limit of massless 6-quarks increase the cross 
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section by about a factor 2 If the QCD corrections to Zjj and Wj] production are of 
similar size, they will weaken the signal to background ratio by about a factor 1.5. 

In the AR and A$ distributions of the system and the jet pair, the Zjj background 
overwhelms the ZZ signal even for rather large ZZV couplings. For the semi-hadronic ZZ 
final states it will therefore be very difficult to utilize these distributions in distinguishing 
the various neutral gauge boson couplings. 

In our discussion of semi-hadronic final states in ZZ production, we have considered 
only inclusive jet rates. With the excellent 6-tagging capabilities of the Tevatron p5| , ^ and 
LHC [^,0 experiments, one may also be able to search for anomalous ZZV couplings using 
ZZ — »• i~^i~bb and ZZ vvhh decays. The ZZ £^£^bb (ZZ —>■ vvbb) cross section is 
about a factor 10 smaller than the ZZ £^£^jj [ZZ i^i^jj) rate if one requires that both 
6-quarks are tagged. The main backgrounds in these channels are Z{^ i~^i~)bb, Z{—^ vv^jb 
and ti production. The Zbb background is about a factor 3 larger than the ZZ i^i^bb 
and ZZ vvbb signals . The signal to background ratio thus is considerably better than 
in the i~^i~jj and z/z/jj cases. However, due to the smaller signal rate, one expects that the 
sensitivity limits obtained from ZZ i~^i~bb (ZZ vvbb ) will be about a factor 2 to 3 
weaker than those derived from ZZ —>■ i~^i~jj (ZZ vvjj\ 



IV. SENSITIVITY LIMITS 



In this section we discuss the method to extract sensitivity bounds on anomalous cou- 
plings, and then determine the bounds on and which one expects to achieve with 
2 fb~^ and 10 fb'^ of data at the Tevatron in Run II, and 10 fb-^ and 100 fb^^ at the LHC. 
For simplicity, we consider only real ZZV couplings. 

We calculate 95% confidence level (CL) limits performing a test. The statistical 
significance is calculated by splitting the selected distribution into a number of bins, each 
with typically more than five events. We use the pT{i~^i~) distribution for all final states 
except vvjj for which we use the prijj) differential cross section. Other distributions, such 
as the ZZ invariant mass distribution (useful only for ZZ 4 leptons), or the maximum 
or minimum transverse momenta of the charged leptons or jets, yield similar results. In 
deriving our sensitivity limits, we combine channels with electrons and muons in the final 
state. In each bin the Poisson statistics is approximated by a Gaussian distribution. Since 
we selected bins containing at least five events, the error introduced by this approximation 
is very small. The same method has been used in the past to estimate limits on anomalous 
WlVj and ZjV couplings for Run I of the Tevatron |1T7| , |42[| . The actual limits obtained 
from experimental data |^-|^ agree well with the predicted bounds. 

In order to derive realistic limits, we allow for a normalization uncertainty of 30% of 
the SM cross section. Backgrounds in the ZZ i^i^vv, ZZ — * C-^d-'jj and ZZ 
vvjj channels are included in our calculation. We impose the cuts described in detail in 



Sec. [Ill A| . In the ZZ t^i~vv case we assume that a jet veto has been imposed to reduce 
the ti background and require pr(^^^~) > 40 GeV to eliminate the Zj background. For 
ZZ jj we require that events do not contain missing transverse momentum of more 

than 20 GeV (40 GeV) at the Tevatron (LHC). Finally, for ZZ vvjj we impose a charged 
lepton veto (pt(^) < 10 GeV if \'r]{€)\ < 2.5) and require exactly two jets to reduce the Wjj 
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and tt backgrounds. As before, we use a form factor of the form of Eq. (j^) with n = 3 and 
A^^ = 750 GeV {ApF = 2 TeV) for the Tevatron (LHC). 

Non-neghgible interference effects between the various ZZV couphngs are found only 
between // and f2, and between and f^. The fY and couphngs do not interfere, 
as expected from the CP-odd nature of /^^ and CP-even nature of f^. This result is 
demonstrated in Fig. [l^, where we show the la and 2a limit contours for ZZ — ^ i'^i'uu 
events in 2 fb~^ of data at the Tevatron Run II. In each graph, only those couplings plotted 
against each other are assumed to be different from their SM values. Plots similar to those 
shown in Fig. |T3| can be obtained for the LHC, other final states, and different values of 
App. As a result of the correlations between and f2, and between and f^, different 
anomalous coupling contributions to the helicity amplitudes may cancel partially, resulting 
in weaker bounds than if only one coupling at a time is allowed to deviate from its SM value. 

In Table | we display 95% confidence level (CL) sensitivity limits expected from the 
Tevatron Run II for integrated luminosities of 2 fb~^ and 10 fb~^, taking into account the 
correlations between different anomalous couplings. Due to the small branching ratio for 
pp ZZ 4 leptons, the number of expected events in that channel for 2 fb~^ is too low 
to allow for an analysis of the Pt{^~^^~) distribution using the method chosen here. The 
bounds obtained from ZZ —>■ i'^i^uu and ZZ vvjj are quite similar. The cross section 
for ZZ — »• vvjj is about a factor 10 larger than that for ZZ —>■ however, the large 

background from Zjj production considerably limits the sensitivity to ZZV couplings for 
ZZ — s> uujj. The limits from the ZZ i^'^^jj and ZZ 4 leptons channels are about 
a factor 1.5 and 2 weaker than those from ZZ l^i'vv and ZZ vvjj. We do not 
attempt to combine limits from different channels. From Table |I| it is clear that this could 
result in a significant improvement of the bounds. 

In Table |T| we display 95% CL sensitivity limits expected from the LHC for integrated 
luminosities of 10 fb~^ and 100 fb^^. The most stringent bounds are obtained from the 
ZZ — > t^l~vv channel. The ZZ —>■ 4 leptons channel yields sensitivity limits which are a 
factor 1.6 to 2 weaker; while, due to the increased Zjj and Wjj backgrounds, the limits 
which can be achieved from analyzing the semi-hadronic channels are a factor 2.5 to 4 worse 
than those found for ZZ —>■ 

The sensitivity bounds on the ZZZ couplings at the Tevatron (LHC) are up to 10% 
(20%) better than those on the corresponding ZZ'y couplings. The sensitivities achievable 
for an integrated luminosity of 10 fb~^ at the Tevatron are about a factor 1.5 to 1.7 more 
stringent than those found with 2 fb~^. At the LHC, increasing the integrated luminosity 
by a factor 10 improves the limits by approximately a factor 2. These numbers are in good 
agreement with the naive scaling law which predicts that sensitivity limits on anomalous 
couplings improve roughly with {JCdtY^^ Due to interference with the SM helicity 

amplitudes, the limits on are somewhat sign dependent. In contrast, because of their 
CP-violating nature, contributions to the helicity amplitudes proportional to fYo do not 
interfere with the SM terms, and the bounds on these couplings do not depend on the sign 
of the coupling. 

The W^W~ and ti backgrounds only marginally influence the sensitivity limits obtained 
in the ZZ t^i~vv case. The limits derived from ZZ ^~^^~jj and ZZ vi/jj 
production, on the other hand, are significantly degraded by the Zjj and Wjj backgrounds 
and thus depend on an accurate knowledge of the signal to background ratio. The signal 
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to background ratio is affected by QCD corrections which are not yet fully known, by the 
factorization scale uncertainty, and by uncertainties originating from the parton distribution 
functions. Changing the signal to background ratio by a factor 1.5 alters the sensitivity limits 
by 10 - 15%. 

The sensitivity limits depend significantly on the scale ApF- For example, at the Teva- 
tron, increasing (decreasing) the form factor scale to App = 1000 GeV {ApF = 500 GeV) 
improves (weakens) the bounds which can be achieved by a factor 1.5 (2). To a lesser de- 
gree, the limits also depend on the power n in the form factor, which we have assumed to 
be n = 3. A smaller value of n allows for additional high px events and therefore leads to a 
somewhat increased sensitivity to the low energy values of the ZZV couplings. 

Limits on the anomalous couplings depend on the power, n, and the scale, App, of the 
form factor. These parameters are a priori unknown, as they represent our ignorance of the 
scale and the nature of new physics beyond the SM. In Ref. it was pointed out that in 
final states without missing transverse momentum one can in principle determine the form 
factor by fitting the distribution simultaneously to f^, App and n. However, a study ||43| 



carried out for WW^y couplings demonstrated that the method will not produce competitive 
limits. One expects that a similar result is obtained for ZZV couplings. Alternatively, 
if non-zero anomalous couplings are observed, the method may be useful in determining 
the shape of the form factor which provides indirect information on the dynamics of the 
underlying new physics. 



V. SUMMARY AND CONCLUSIONS 

ZZ production in hadronic collisions provides an opportunity to probe neutral gauge 
boson self-interactions in a direct way. In this paper we presented a detailed investigation 
of how well future experiments at the Tevatron and LHC will be able to measure the ZZV 
couplings. Our calculation has been carried out in the double pole approximation, using 
the most general ZZV vertex function which respects Lorentz and electromagnetic gauge 
invariance. If both Z bosons are on-shell, there are two independent ZZZ (// and f^) 
and two ZZj couplings {f2 and f^). Z decays with full decay correlations and finite Z 
width effects are included in our calculation. Non-resonant Feynman diagrams, except for 
time-like photon exchange diagrams, and contributions from gg —>■ ZZ are not taken into 
account. The non-resonant Feynman diagrams change the differential cross section by at 
most 5%. The contribution from gluon fusion enhances the cross section by about 1% (15%) 
at the Tevatron (LHC). NLL QCD corrections are approximated in our calculation where 
appropriate by a simple K-factor. 

S'-matrix unitarity requires that the ZZV couplings are momentum dependent form 
factors. We derived constraints on the low energy values of the ZZV form factors and 
the shape of the form factor from a partial wave analysis of the inelastic ZZ production 
amplitude in fermion antifermion annihilation (see Eqs. ( p!0D and ([TT|) ). 

The effects of anomalous ZZV couplings are enhanced at large energies. They lead to a 
broad increase in the ZZ invariant mass distribution and in various transverse momentum 
distributions. We considered four signatures of ZZ production: ZZ —>■ 4 leptons, ZZ — »• 
i~^£~ui', ZZ —>■ and ZZ vvjj. The ZZ —>■ 4 leptons channel is almost background 

free but suffers from a small event rate due to the small branching ratio for a Z pair decaying 
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into four charged leptons. The rate for ZZ t^i~vv is approximately a factor six larger 
than that for ZZ 4 leptons; however, W~^W~ and ti production constitute non-negligible 
backgrounds at small values of pT{i~^i~)- The ZZ C-^C-' jj and ZZ vvjj modes have 
larger branching ratios than the ZZ 4 leptons and ZZ t^t~vv channels, but also 
much higher backgrounds. The main background source for the semi-hadronic final states 
is QCD Z + 2 jet production. The Z + 2 jet cross section is about a factor 20 (50) larger 
than the ZZ signal at the Tevatron (LHC). Nevertheless, we found that these channels may 
be useful in searching for ZZV couplings. 

All ZZV couplings have similar high energy behavior. This makes it difficult to dis- 
tinguish the various ZZV couplings in the transverse momentum and ZZ invariant mass 
distributions. We found that in the ZZ 4 leptons and ZZ t^i~vv channels, the 
distributions of the lego plot separation, l\R{Jt^l~\ and the angle in the transverse plane 
between the charged leptons, A$(£"'"£~), may be useful to distinguish fX from and to 
determine the sign of Since violate CP conservation, terms in the helicity amplitudes 
proportional to do not interfere with the SM helicity amplitudes, thus it is impossible 
to determine the sign of fX using the Ai?(^"'"^~) and A$(£+£~) distributions. In the semi- 
hadronic modes, the ZZ signal is overwhelmed by the Zjj background in these distributions, 
even for rather large anomalous couplings. Unfortunately, the two distributions will not be 
useful in distinguishing between ff and as they differ very little in shape for ZZZ and 
ZZ^I couplings. 

In order to determine the bounds on the ZZV couplings which one can hope to achieve in 
future Tevatron and LHC experiments, we have performed a test for the four different final 
states, using the t^t~ and dijet transverse momentum distributions. At the Tevatron, with 
an integrated luminosity of 2 fb~^, one will be able to measure j\ and with a precision 
of 15 — 20% (95% CL) in the ZZ — > i'^i^vv and ZZ vvjj channels for a form factor 
scale oi KpF = 750 GeV. The limits obtained from the ZZ 4 leptons and ZZ ^^^^jj 
channels are a about a factor 2 and 1.5 weaker than those which can be achieved in the other 
two channels (see Table |). At the LHC with 100 fb~^, the ZZV couplings can be measured 
with an accuracy of 3 x 10^^ - 4 x 10^^ (95% CL) in pp^ ZZ ^ for A^.^ = 2 TeV. 

The limits obtained from the other three channels are a factor 1.6 to 4 weaker (see Table y). 

The sensitivity limits which can be achieved at the Tevatron in Run II and the LHC 
should be compared with the bounds from recent measurements at LEP2 and expectations 
for a future linear collider, as well as predictions from theory. The combined 95% CL limits 
from the LEP2 experiments presently are |^ 



- 0.66 < /f < 0.68, -0.40 < fj < 0.38, (17) 
-1.06 < < 0.69, -0.89 < f2 < 0.86. (18) 

It should be noted that the LEP2 limits do not contain any form factor effects. For the 
form and scale chosen here, these effects weaken the limits by about 20%. In Run II, CDF 
and D0 will be able to improve these bounds by at least a factor 3 to 6. At a e^e~ Linear 
Collider (LC) with ^/s = 500 GeV and an integrated luminosity of JCdt = 100 fb~\ one 



expects to measure the ZZV couplings with a precision of |45| 4x10 — 6x10 (95% 



CL). This is comparable to the bounds which one hopes to achieve at the LHC. One loop 
corrections in the SM, in supersymmetric models, and in models with heavy fermions, induce 
ZZV couplings which are of O{10~^) or less. 
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In view of our present poor knowledge of the ZZV couplings, the direct measurement of 
these couplings in Run II and at the LHC will constitute major progress. However, it will 
be very difficult to achieve a precision which will make it possible to test the SM one loop 
prediction for the ZZV couplings. 
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TABLES 



TABLE L Sensitivities achievable at 95% CL for anomalous ZZV couplings in 
ZZ 4 leptons, pp ZZ l^i^uu, pp ZZ ^^^^jj, and pp ZZ vujj at 
the Tevatron {y/s = 2 TeV) for an integrated luminosity of (a) 2 fb~^, and (b) 10 fb~^. The limits 
for each coupling apply for arbitrary values of the other couplings. For the form factor we use the 
form of Eq. (|2|) with n = 3 and KpF = 750 GeV. The cuts imposed are discussed in the text. 

(a) jCdt = 2 fb-i 



coupling 


ZZ 


4 leptons 


ZZ i+l'uv 


ZZ i+t 


33 


ZZ vvjj 


/40 






+0.169 
-0.169 


+0.219 
-0.218 




+0.159 
-0.160 


/40 




— 


+0.175 
-0.174 


+0.222 
-0.221 




+0.163 
-0.162 


/so 






+0.171 
-0.204 


+0.220 
-0.244 




+0.162 
-0.184 


/50 






+0.184 
-0.202 


+0.229 
-0.241 




+0.170 
-0.179 


coupling 


ZZ 


4 leptons 


(b) ^Cdt = 10 fb-i 
ZZ i+l-uu 


ZZ e+t 


'33 


ZZ vvjj 


/40 




+0.180 
-0.179 


+0.097 
-0.097 


+0.146 
-0.145 




+0.106 
-0.106 


fl 




+0.185 

-0.185 


+0.100 
-0.099 


+0.148 
-0.147 




+0.109 
-0.108 


/so 




+0.178 
-0.216 


+0.092 
-0.120 


+0.144 
-0.167 




+0.105 
-0.127 


/so 




+0.192 
-0.213 


+0.103 
-0.115 


+0.151 
-0.163 




+0.112 
-0.121 
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TABLE II. Sensitivities achievable at 95% CL for anomalous ZZV couplings in 
•pp — > ZZ — > 4 leptons, — > ZZ pp ZZ and pp — > ZZ — > v'ujj at 

the LHC {y/s = 14 TeV) for an integrated luminosity of (a) 10 fb"\ and (b) 100 fb"^ The hmits 
for each coupling apply for arbitrary values of the other couplings. For the form factor we use the 
form of Eq. (|2|) with n = 3 and App = 2 TeV. The cuts imposed are discussed in the text. 

(a) jCdt = 10 fb-i 



coupling 


ZZ ^ 4 leptons 


ZZ i+l'vu 


ZZ ^ l+l-jj 


ZZ uujj 




+0.0115 
-0.0114 


+0.0060 
-0.0060 


1 n noon 

+0.0230 
-0.0228 


+0.0156 
-0.0154 




1 n m on 

+o.oioy 

-0.0139 


+0.00 / 2 

-0.0072 


+0.02/4 
-0.0274 


+0.0186 
-0.0186 


/50 


+0.0119 
-0.0113 


+0.0052 
-0.0060 


+0.0226 
-0.0220 


+0.0160 
-0.0158 


fj 

J 50 


+0.0137 
-0.0145 


+0.0072 
-0.0075 


+0.0268 
-0.0274 


1 n m oo 

+0.0188 
-0.0190 






(b) JCdt = 100 fb-i 






coupling 


ZZ 4 leptons 


ZZ l+l-uv 




ZZ ^ Pzvjj 




+0.0052 
-0.0051 


+0.0031 
-0.0031 


+0.0130 
-0.0128 


+0.0088 
-0.0086 


fy 


+0.0062 
-0.0062 


+0.0038 
-0.0038 


+0.0154 
-0.0154 


+0.0104 
-0.0104 


/so 


+0.0053 
-0.0051 


+0.0032 
-0.0031 


+0.0128 
-0.0122 


+0.0092 
-0.0088 


/so 


+0.0061 
-0.0065 


+0.0037 
-0.0039 


+0.0148 
-0.0154 


+0.0104 
-0.0108 
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FIGURES 




FIG. 1. The Feynman diagrams for the tree level processes contributing to pp ZZ in the SM. 




q Z 



FIG. 2. Contributions of ZZZ and ZZ^ diagrams to qq ZZ. 




FIG. 3. Feynman rule for the general ZZV iy = Z, 7) vertex. The vertex function is given in 
Eq. (H) . e is the charge of the proton. 
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FIG. 4. The pt{Z) and PtIm) distributions in pp ZZ — > e~^e /i^/i at the Tevatron 
{^/s = 2 TeV) in the SM (sohd Une), for /J) = 0.3, and f^Q = —0.3. The cuts imposed are 



described in detail in Sec. Ill A. The form factor scale has been set to App = 750 GeV. 



25 



10^ 



> 
o 



10" 



10" 



10" 



i 10-4 



b 
t3 



10" 



10" 



10" 



I" " 1 ' 

: a) 


" 1 " " 1 " " 1 " " 
pp ^ ZZ ^ e~ jj,'^ jj,' 




- \ 


Vs = 14 TeV 


- 


- w 








N 

- . \ 
" - \ 
■ .\ 

\ ^ 
\ ■ ^ 
\ • ^ 
\ ^ 




: solid: 
dash: 


\ ■ \ 
\ ' ■ 

SM \ 
f4o = 0.02 


\ 

■ ,\ 

\ = 
■ ■ .\ I 


: dots: 

, , , 1 , 


f?o = -0.02 

, , 1 , , , , 1 , , , , 1 , , , , 


1 



10^ 



10" 



10" 



10" 



10" 



10" 



250 500 750 1000 1250 1500 
Pt(Z) (GeV) 



10" 



b) pp^ZZ^e^e 
Vs = 14 TeV 




200 400 600 800 

Pt(m) (GeV) 



FIG. 5. The pt{Z) and pxifJ-) distributions in pp ZZ e+e /i^/i at the LHC 
= 14 TeV) in the SM (solid line), for /|, = 0.02, and f^^^ = -0.02. The cuts imposed 



are described in detail in Sec. [II A. The form factor scale has been set to App = 2 TeV. 
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FIG. 6. The AR{ij,^fi ) and A^{ij,^ ) distributions in pp — > ZZ e^e ^^/.i at the Teva- 
tron (y/s = 2 TeV) in the SM and in the presence of non-standard ZZZ couphngs. The cuts 



imposed are described in detail in Sec. Ill A. The form factor scale has been set to App = 750 GeV. 
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FIG. 7. The AR{fj,^ ) and A^^n'^fi ) distributions in ZZ e^e fi^fi at the LHC 

{^/s = 14 TeV) in the SM and in the presence of non-standard ZZZ coupUngs. The cuts imposed 



are described in detail in Sec. [II A. The form factor scale has been set to App = 2 TeV. 



28 



10^ 



> 



E-1 

^ 10-2 
b 



10" 



"T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



a) 



pp ^ e e w 
Vs = 2 TeV 




dash: SM full 

solid: SM ZZ + W+W" appr. 



1.15 



1.10 



^ 1.05 

a 

ft 

(0 

1.00 



0.95 



0.90 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



b) pp^e'''e uu 
Vs = 2 TeV 



0.85 

50 100 150 200 50 100 150 200 250 300 



PT(e+e ) (GeV) 



PT(e + e ) (GeV) 



FIG. 8. Full and approximate results for the transverse momentum distribution of the e+e~ pair 
in pp e^e~uu at the Tevatron {^/s = 2 TcV) in the SM. The individual distributions obtained 
using the full set of contributing Feynman diagrams (solid line) and the subset of diagrams for 
qq ZZ — > e^e'vu and qq —>■ W^W~ e'^e'Uf.i'e in the double pole approximation (dashed 
line) are shown in part a). Part b) displays the ratio of the full and approximate differential cross 
sections. The cuts imposed are described in detail in the text. 
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FIG. 9. Transverse momentum distribution of the e+e~ pair in pp ZZ e^e'vu at the 
Tevatron (-y/i = 2 TeV) for the SM and for = 0.3, together with the differential cross sections 
from several background processes (a) without and (b) with a jet veto applied. The cuts imposed 
are described in detail in the text. The form factor scale for nonzero ZZV couplings has been set 
to Aff = 750 GeV. 
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FIG. 10. Transverse momentum distribution of the e'^e~ pair in pp ZZ e^e~vu at the 
LHC ("v/i = 14 TeV) for the SM and for = 0.02, together with the differential cross sections 
from several background processes. The cuts imposed are described in detail in the text. The form 
factor scale for nonzero ZZV couplings has been set to AiT'p = 2 TeV. 
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FIG. 11. The Ai2(e+e") and A$(e+e") distributions in pp — > ZZ e^e vv at the Tevatron 
{^\fs = 2 TeV) in the SM and in the presence of non-standard ZZZ couphngs. The cuts described 
in Sec. Ill A and a > 20 GeV cut are imposed. In addition, we require that no jets with 
PtO) > 20 GeV and \'r]{j)\ < 3.5 are present. The form factor scale has been set to App = 750 GeV. 
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FIG. 12. The Ai?(e+e~) and A^>(e+e~) distributions in ^ ZZ ^ e^eTvv at the LHC 
= 14 TeV) in the SM and for /I, = 0.02 with a form factor scale of A^^i^ = 2 TeV. The dotted 
and dash-dotted curves represent the tt and W~^W~ + jet backgrounds. The cuts described 



in Sec. Ill A and a ;^rp > 50 GeV cut are imposed. In addition, we require that no jets with 
PtU) > 50 GeV and \r]{j)\ < 5 are present. 
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FIG. 13. Transverse momentum distribution (a) of the £'^£~ pair in pp — ZZ — £'^£~jj, and 
(b) of the jet pair in pp — >■ ZZ vvjj, at the Tevatron {^/s = 2 TeV). The SM prediction is shown 
together with the cross section for = 0.3 and Kpp = 750 GeV. Also shown are the differential 
cross sections of various background processes. The cuts imposed are described in detail in the 
text. 
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FIG. 14. Transverse momentum distribution (a) of the £'^£~ pair in pp — ZZ — £'^£~jj, and 
(b) of the jet pair in pp — > ZZ i^i^jj, at the LHC {^/s = 14 TeV). The SM prediction is shown 
together with the cross section for = 0.02 and App = 2 TcV. Also shown are the differential 
cross sections of various background processes. The cuts imposed are described in detail in the 
text. 
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FIG. 15. Correlated sensitivity limits for ZZV anomalous couplings in pp ZZ t^fvv 
events at the Tevatron Run II (-y/s = 2 TeV) with 2 fb~^ of data. All couplings are assumed to 
be real. Shown are la (dashed lines) and 1g (solid lines) limit contours of all combinations of jf^ 
versus /^^ (i = 4,5). A form factor scale of h.pp = 750 GeV has been assumed. In each graph, 
only those couplings which are plotted against each other are assumed to be different from their 
zero SM values. 
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